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M  itli  <|H>>4Mt'«»4 


a'.w  ^  fNMp  «'|lt«il»'  imP  ft JNifi«liirl*»  •»!  IH*# 

.^.a  «w.a>ai4'4  iMk^lliMf*  git  '4Nr|g«Mg>  iMal  ftMNMlIC^  Hmm 

IMI#  iHig  w«ri«NHi  ii«ff€i«  tut 

iit<u«i  «t.'0up  •I«i4l«»  giw«g  l«%7.  Hilg  t«MMr<fi  Imm 

jixH^  8Wii«»#||<f4l  gl«g4l|g«»  figlll  (MNMMIVgllMI  |»r««rfllt«  iAia  waif* 

*r«,  «m4  iibi«o  4«^igp»  «i«|  gm«lriK'il«n  «f  apaelal  a^alpMai  lo  laflanMl 
tw  (erl<l  «t4  l4i«*rgt«rv  fiM|ttlrgMMta* 

p«f«  iiM^rift  iMg  MM  Mfongorvil  la  tlia  ^aac  kjr  tha  tl>*  R.  Aragr  Rlipuil 
«««rfig,  tSw  («rfi€«  of  Xaval  Raggareh*  tha  Air  fbrca  CanArl^ita  Raaaarch 
luthitrotorlf^,  tha  Air  fbrca  Tachnlcal  il^llcaKlaaa  CaaKar»  cha  Air 
I'orctf  offitrg  of  Sclantlflc  Raaaarch  and  cha  NatlMal  Relaaea  fbuada- 
tion«  Tha  work  raportad  haraln  covara  a  four  yaar  klM  pariod  froa 
1  Juno  Hftit  elirrufth  31  May  1970  and  waa  aupportad  by  cha  Advancad 
Ranonrcl)  Projccca  Aaancy,  ARPA  Ordar  No.  292,  AMndMntt  32  and  37, 
and  W4ifi  non!  corad  by  cha  Air  Porca  Of  flea  of  Sclonciflc  Raaaarch  undar 
C-ontract  AF4»»(63R)-1759, 

Thin  report  la  wriccen  Co  fulfill  conCracCual  raqulrananCa  for  a 
final  report.  Previous  technical  reports  and  scientific  Journal  arCi> 
cles  have  presented  a  large  portion  of  the  research  sponsored  by  this 
contract,  and  the  data  contained  in  these  reports  and  articles  will 
only  be  referred  to  In  this  report. 
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•lli»>«l%«W(°f 

4^  fftoW  ’4MF|A-Iji»|l  mt  «SI4tH*S.<*l  «iKS«4li4  *  fi*  «AKM 

.‘ifVfWIltli^Nil  l#)f  ‘«)lillll»¥4lM  Sit4  4M|(' @«»y«<^'lo 

llli|ll4i  wVMIpiNSVf  I*  tm  -4  «4^#g9i-«i9  mi  4  t$4l  •«44il4'<0 

(if  t>t(f44-»«l>4IW*kK(iiill  (Wxiplat  4f4  4M  #4>4'||<)>H4»  wf  |l<t-»  4#» 

«lifii««l  tiif  |i(t  wlfW'llf  illMfl^iilliMi  l««  ^14  #c|»»»<*«)4i4v«  4g|4««  4*  m 
4Mlffl(r«l  flHKfetiMl  «f  C«I  «  tflil.  ii^4|'4  4f4  firiivtl  !Hl  OWpop 

l«f  Mil  iniwl  IliMMli  4f«  C4l(r«ll4l#4l  tc»  Mf  tlifilftpd  «l»ftr44  tif  4r«wf4f'f 
for  IlM  Motfol  mdif  (mosl^orotloo.  Hi#  iMrllMMi  coo  to#  mtinf  tm  thmirnfm 
teoioow  Mo^U,  Mul  laeloiM  towpoiotioo  of  nolcl^l#  roflociloM. 
Trovol  Kioto  hoot  htoo  ohiolood  for  o  nodtl  of  o  (ffooo-oociloo  ootltf- 
lying  LASA«  In  o  auyyltKtae  (flOf|*))*f|)  lo  chin  finol  r»|M>re  Ihin 
oiMiloelon  ntehod  la  donerlbod  and  llluoiroctd  In  dtcoll.  Hm  nupyl#- 
ncne  la  btaod  upon  a  doecorni  dlaoortaclon  on  thin  ntchod  hy  f.  U 
Jaekaon* 

In  1467  cht  Invta cl nation  of  opcieal  proctaolnt  waa  canalnacadt 
tniciacad  In  1461,  cht  Inotaclgatioo  oaa  conducctd  Co  txplort  cht 
appllcaClon  of  coharanc  opcical  procttaing  Co  atlanlc  analyala.  Ica 
uat  waa  found  Co  lit  primarily  in  oil  txploraclon  atianoloRV*  Tli# 
reaulca  of  cht  opClcal  invtaclRaCion  wtrt  proatnCed  In  atvtral  journal 
arCiclea  during  Cht  period  1963-1968. 

A  daca  proctaaing  capability  apecifically  inCtnded  for  the  atudy 
of  Che  S-waye  parCiclt  variaCiona  acroaa  LASA  waa  built  on  thia  con¬ 
tract.  The  LASA  data  processor,  which  existed  in  this  laboratory  for 
the  optical  analyses  of  LASA  data  was  modified  to  interface  with  this 
laboratory's  PPP-8  computer.  This  would  allow  unlimited  computinR  at 
no  cost  and  inmediate  man-machine  interaction.  An  inhouse  wiggly  line 
output  device  was  modified  to  print  out  the  processed  and  unprocessed 
three-component  data.  This  capability  was  not  used  due  to  funding 
limitations. 

Microseismic  investigations  in  connection  with  JORUM  indicated 
that  the  firing  of  a  large  yield  underground  shot  did  not  produce  any 
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•*i'  tll'H’  t**  '  H  fJItl  ,.(  a  MV  P 

iv,, .vj  M  4.  -jiii  4^  iiiA.,Ht  «.!i»  -♦w'i!n#M!i,» 

.  H^tlll  Ite04«*44*fltt4l4>  4v««>  1*«#  Srt#'  iflta 

-t,  *-»4<4.(,M'' =  01  .a  ^'ki  ^h**#*'  fs# 

iv.«  4lffi*<|«o  ^  IiIh 

.(•kacJHN,  4IHf  ^$»m-  9(t^»  «Ni»  mwimil 

ffliltJ-  #ik^4>l  4Fk«^  $4  4I¥i  9(t4  #<IWH>il><^»i  I<lilft4l| 

.t4«i.<lj,»444g,4  nig  S.'4>H»e 

Hat4w«-H  •»4  4«tll#>*<  x'-'f  .'•H#  «4Pf«>  iMtl  ••  dMftvimM*  tA#  *««#« 

(  4uw>  8«H4  #if  «  fiM»f1|i(MI  IM»  MiftA)*#  •<••••  All  MIIMMilW 

«'.«*  itt««f>  m  fi»  4ii«>Miiil(Mi  mi  ttmml  liMm  at  Mlfliilc  mmm 

IM  llo««l##W  HDtllctl  t1lM»f»l  l»4ll  ilwilM  Wtf*  mA»  ••  4v» 

(••tmiiffir  iM  •!  fl«i9  •latilc  VMM  liy  attrfM* 

|irrfv<fili««*  i«M  i««  «ll»i#mtii#  iImp  mt  •rr*r  l«  ilit  •l••-l•nl 

In  nnlMfMlOKf* 

Spnrirnt  «fi»(ll«<>«  wvr*  n*ln  of  WNlofiroMoA  Boeloor  bImco  ami 
»4«rciit«|M.i)(»9i  r«*r<iriln«l  at  l•l•••l•llic•  roiloool  ood  otar  diotic— » 

An  <i«l<llci««n<il  niudy  won  niidi  of  Clio  clinraecorioclcs  of  ooImIc  bock* 
i:rirtun«l  nolnn  ro^ordnd  on  cho  HocCm  of  t.oko  Siiportor* 
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Hm  D  %MtN4l^MN«  n>«t  «Mj 


IpM  llRp»*«lfp4l  Itt  414^141(11  i«Mv 

(lUm  *(I^  I#  |A#9#  lW'l««llf  I4iv«il«fv4  41  fftt# 

^4»44Hlt4  IWMiV  4m  |4«)  |k»  irAt»l|«fV  »t 

4  im  in*#  fvIvfiWifv  UtlM'Mf  «iMh|»K«  wm  mrnU  l«9  llw 

It4i«ll4l  IIUNP  At  Il»  |llif|Pl4llR«l  f!«»‘  mt  11^4  4«4Vi«lFf<iit<W  |i«  ^  ttm>  ImrlilMI 
.4inAlA  At  PIjAA  A4iA  I A  IftA  *  MhtfAA  |4«#*.  tW  4«4(|^4lA»  ifiMII  AM  I  .I 
<fWMA  fAf  (NKA(  fmif  AHUPMA^c  l«MrtiMllA«  tflAI  At4  mf  IaaMIam  At  9m  f 
mm<0>  At  fiiqp,  CAMVl  IldA,  4M  APPfAAlnAtA  fAl4li«A 

oimrM  AMLVSis  tf«iiii<iti.« 

A  tilfti#)  latNMtliMtlM  At  trA(«Mi««  fAr  mIamIc  iaia 

mm  «8ac.liii»i  la  lAAl*  Hila  lavaatliailaa  aaa  laltlaia^  la  lAAl  aaaar 
AHWA  fiMifatt  Af  tmi  caaelaaatf  aaAar  AIQAi  CDacracK 

AH<KAia|»iilA«  Hiia  laatatlAailaa  aaa  ilia  firai  la  i^lf  caiiaraat 
apilcai  prataaalai  la  aalaaica.  It  8ilaalaia4  mmf  athar  la«aatlMilMfi 
alililli  fMaica4  la  th9  aptlleatlaa  af  aftlcal  ^laeaaalaa  to  all  osfilora- 
tiaa  aaliaalny,  md  ciw  MflwclaA  of  caaaorclal  laatnanwta  for  audi 
fracaailai*  Haaaaar,  lia  iptllcailoa  to  VILA  tMIftNOt  ^roAloat  aao  not 
rafflcloatlf  aMfal  ta  caatlaat  tlw  lavoatliatlaa.  Ofitlcal  ^rocaaolmi 
aaa  aaAt  loaa  attractloa  bjr  tlM  atfvaat  of  tko  foot  fburlor  tranofoia 
prasrana  far  tlio  digital  coagator* 

At  tho  Incagtioa  of  thla  contract  In  l%A  it  aaa  plwnod  to  apply 
Iwth  eaMfcnt  and  non-cohomat  optical  procanalng  to  Montana  lASA  data* 

A  data  convertor  waa  conatrvetod  to  aoloct  data  channolo  fron  LASA 
digital  tapoa  for  recording  on  varlablo-danvity  filna.  Ihia  convertor 
wao  oonotmetod  and  filaa  aere  nado*  However*  tho  deciaion  by  the 
Sponnor  to  create  tho  Seiaalc  Analyai*  Cantor  nado  tho  optical  approach 
ohaoloto  In  tho  opinion  of  tho  principal  invoatigator  and  the  Sponaor. 
lha  inveatlgation  wao  therefore  teminatod.  Tho  LASA  data  convertor 
haa  bean  uaad  aubooquantly  to  aoloct  and  ro-rocord  I  ASA  data  for  other 
Invoati gat Iona  within  tho  Goophyoica  Laboratory. 
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uliiii*  till*  •'lUii r.ict  |)t•rlOil  wt«rk  Iuim  prottruaHed  toward  an  under* 

At .isi.liin;  Ki  *^-v.nv  particUi  notion  variattona  acroas  LASA*  Due  to 
li^vitoii  fniKln  li  wiin  iloclileil  to  do  the  nocoaaary  aignal  analysis  on 
till'  dtitit.il  conputor  that  is  in  our  laboratory  and  available 

.It  .ill  tini'fi  at  no  cost.  Altliougli  this  conputer  is  much  smaller  than 
till'  litM  ii<n/K7  available  on  campus  it  offers  the  additional  advantage 
(ntiii'i  than  no  cost  I  of  irnncdiatc  man-machine  Interaction  even  when 
the  data  is  >>n  tape. 

The  input  for  the  IMH'-H/I  is  the  l^A  data  processor  and  it  had 
to  be  interfaced  with  the  I’DP-A/I*  Tliia  work  has  proceeded  wall  and 
is  'lescrihed  below. 

An  interface  was  designed  and  constructed  between  the  LASA  data 
converter  and  the  PDP-H/I  digital  computer.  Via  this  interface,  under 
computer  program  control.  LASA  digital  data  can  be  stored  in  the  PDP- 
H/i  tiemory  for  digital  processing.  The  LASA  data  converter  assembles 
the  digital  seismometer  samples  from  the  bytes  on  magnetic  tape  and. 
at  the  proper  time,  the  computer  accepts  that  sample.  Since  the 
IMM'-S/l  IS  a  12-hlt  word  length  machine,  a  switch  setting  selects  the 
12  most  significant  bits,  the  12  least  significant  bits  or  the  "middle" 
12  bits  of  the  14-bit  LASA  word. 

Mata  from  either  LASA  high  rate  or  low  rate  tapes  may  be  stored 
via  the  interface.  A  maximum  of  SO  words,  in  either  format,  may  be 
selected  on  the  IJ\SA  data  converter  and  those  words,  taken  in  increas¬ 
ing  numerical  ord^r  in  each  record,  will  he  stored  in  sequential  ad¬ 
dresses  in  the  PMI’-8/I  memory.  The  number  of  samples  to  be  stored 
may  be  controlled  either  by  the  conq>uter  program  or  by  the  number  of 
records  the  LASA  data  converter  reads,  whichever  method  is  appropriate. 
Available  computer  memory  space  governs  the  total  number  of  LASA  sam¬ 
ples  that  may  be  stored. 
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To  check  the  total  system  LASA  data  tapes  were  required.  Since 
great  delays  were  encountered  In  receiving  these  tapes  as  has  been 
discussed,  a  request  was  made  directly  to  M.  Nafi  Toksoz  of  Planetary 
Sciences  Department  at  M.I.T, ,  to  get  a  composite  tape  and  the  S  waves 
used  by  John  W,  Fairborn  (1969)  in  recent  paper.  This  data  was  re¬ 
ceived  promptly  and  the  total  LASA  data  processor/PDP-8/I  system  is 
now  ready. 

No  output  device  displaying  wlggly-llnes  was  available  for  the 
PDP-A/I.  Again  due  to  limited  funds  none  could  be  purchased  on  this 
contract,  therefore,  Inhouse  funds  are  providing  a  three-channel  ana¬ 
log  output  capability  for  the  system. 

Despite  the  fact  that  this  data  processing  system  Is  now  ready 
for  use  at  no  cost,  work  was  terminated  due  to  depleted  contract  funds. 
Additional  funding  would  be  required  to  complete  the  anticipated  work. 
Since  we  are  only  now  ready  to  begin  the  actual  S-wave  analysis  and 
since  the  necessary  computing  can  now  be  done  at  no  cost  It  Is  sin¬ 
cerely  hoped  that  additional  funding  can  be  obtained. 

2.4.  SPECIAL  STUDIES 

2.4.1.  MICROEARTHQUAKE  INVESTIGATIONS.  The  Geophysics  Laboratory 
operated  a  tripartite  net  of  three-component  short-period  seismograph 
stations  for  a  period  of  two  weeks  before  the  JORUM  underground  nuclear 
shot  and  for  three  weeks  after  to  study  the  mlcroselsmlc  activity  In 

a  tectonically  active  area  northeast  of  the  shot.  The  results  of  this 
study  are  contained  In  a  paper  entitled,  "Mlcroearthquake  Studies  In 
Connection  with  the  JORUM  Underground  Nuclear  Detonation,"  by  D.  E. 
Willis  that  has  been  submitted  to  the  Bulletin  of  the  Selsmologlcal 
Society  of  America  for  publication. 

2.4.2.  DECOUPLING  EXPERIMENT.  Five  short-period  seismograph  sta¬ 
tions  were  operated  at  the  Nevada  Test  Site  to  record  DIAMOND  DUST,  a 
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low  yield  nuclear  device  that  was  detonated  in  a  spherical  cavity  to 
determine  the  feasibility  of  using  thermomechanical  techniques  for 
decoupling  seismic  energy.  Excellent  results  were  obtained  from  this 
experiment.  A  technical  report  describing  these  results  will  be 
published  as  a  supplement  to  this  report  (8071-33-F3/3620-1-T)  that 
will  present  the  results  of  this  test. 

2.4.3.  CONFIGURATION  OF  THE  SURFACE  OF  THE  EARTH'S  CORE  FROM 
PcP  TRAVEL  TIMES.  A  comprehensive  study  (doctoral  thesis— C.  Bufe) 
was  made  using  PcP  and  P  travel  times  to  ascertain  the  configuration 
of  the  surface  of  the  earth's  core.  Nuclear  and  high  explosive 
sources  were  used  in  this  study.  Travel-time  residuals  from  the 
Taggart-Engdahl  (1968)  tables  were  determined  for  PcP-P  time  Intervals 
and  for  PcP  travel  times,  corrected  for  station  and  source  terms. 
Corrections  were  also  made  for  elevation  of  source  and  station  and 
for  ellipticity,  excluding  the  core  ellipticity  term. 

The  distribution  of  the  data  would  not  allow  the  application  of 
conventional  spherical  harmonic  analysis  techniques  on  the  Inferred 
variation  of  the  core  radius.  Hence  a  modified  spherical  harmonic 
analysis  method  was  devised  to  smooth  the  data  and  to  estimate  the 
shape  of  the  core.  The  variation  of  the  core  radius  was  determined 
to  be  ±10  km  using  this  technique. 

This  investigation  is  described  in  detail  in  a  supplement  (8071- 
33-F2)  to  this  final  report.  The  supplement  is  based  upon  a  doctoral 
dissertation  by  C.  G.  Bufe. 

2.4.4.  OTHER  STUDIES.  A  number  of  other  studies  were  conducted 
during  the  course  of  this  contract.  These  have  been  reported  in 
journal  articles,  University  of  Michigan  technical  reports  and  at 
scientific  meetings.  These  reports  are  listed  below  and  abstracts  are 


6 


TVUt^ 


iNSTirtiTi';  or  sciinci:  and  riciiNoixuiv 
nir  I  Nivi  Rsn  v  or  miciik.an 


reproduced  in  Appendices  I-III : 

"A  Long  Range  Reversed  Seismic  Refraction  Profile 
in  the  Eastern  U.  S."  by  D.  E.  Willis  and  J.  T.  Wilson. 

"Seismic  Radiation  from  SHOAL"  by  C.  G.  Bufe  and 
D.  E.  Willis. 

"Seismic  Noise  on  the  Bottom  of  Lake  Superior" 
by  D.  E.  Willis. 

"An  Investigation  of  Seismic  Wave  Propagation  in 
the  Eastern  U.  S."  by  D.  E.  Willis. 

"Scattering  of  Plane  Elastic  Waves  by  Surface 
Imperfections"  by  I.  K.  Mclvor. 

"High  Frequency  Teleseismic  Energy  From  Aleutian 
Sources"  by  C.  G.  Bufe  and  D.  E.  Willis. 

"An  Investigation  Into  the  Time-Term  Method  in 
Refraction  Seismology"  by  Leon  Reiter. 

"A  Note  on  Short-Period  Seismic  Spectral  Ratios" 
by  C.  G.  Bufe  and  D.  E.  Willis. 

"Near  Earthquake  Spectra"  by  C.  G.  Bufe  and 
D.  E.  Willis. 

"A  Note  on  the  Anna,  Ohio  Earthquake  of  July  26, 
1968"  by  D.  E.  Willis  and  James  T.  Wilson. 


3 

CONCLUSIONS  AND  RECOMMENDATIONS 
3.1.  DIGITAL  RAY  TRACING 

A  versatile  tool  has  been  developed  for  the  scientific  investiga¬ 
tion  of  velocity  distributions  in  the  earth.  Any  conceivable  velocity 
distribution  can  be  modeled,  and  by  simulation  the  travel  times  can 
be  compared  with  actual  travel  times  to  determine  the  validity  of  the 
model.  The  simulation  applies  wherever  geometrical  ray  approximations 
are  useful. 
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It  is  ri'commendecl  that 

1  .  The  programs  he  made  as  usable  as  possible  for  other 
investigators  and  funds  for  providing  cards,  programming  help,  and 
dissemination  to  other  Investigators  be  provided. 

2.  An  extensive  investigation  of  models  of  the  underlying 
structure  of  LASA  and  other  pertinent  regions  be  Initiated. 

3.  The  method  be  extended  to  Include  mode  conversion  and 
/'.oepprltz's  equations. 

4.  Travel  times  and  positions  are  known,  therefore  the 
basic  information  for  computing  dlff  "Ive  effects  is  available,  and 
should  be  investigated, 

3.2.  OPTICAL  ANALYSIS  TEOINIQUES 

This  investigation  has  been  terminated. 

It  is  recommended  that  no  further  effort  be  made  in  this  direction. 
Subsequent  to  the  investigation  reported  here,  there  have  been  many 
investigations  of  optical  analysis,  primarily  by  the  oil  exploration 
companies.  Duplicating  these  efforts  would  be  redundant.  One  exten¬ 
sion  of  this  technique  is  holographic  seismology,  which  has  been  in¬ 
vestigated  in  three  laboratories  for  approximately  five  years.  Results 
will  be  disseminated  if  holography  is  useful  for  seismics.  In  the 
opinion  of  the  principal  investigator  there  is  little  hope  that  holo¬ 
graphy  will  be  useful,  except  possibly  heuristically  in  giving  new  in¬ 
sights  which  might  result  in  new  digital  computation  techniques, 

3.3.  MODE  FILTERING - S-WAVE  STUDIES 

Although  several  array  studies  (Fairborn,  1969;  Whitcomb,  1969)  of 
S  waves  have  been  made  it  appears  that  no  work  has  been  published  con¬ 
cerning  explosion  generated  SH  waves  propagating  across  LASA,  Work 
currently  being  done  by  Turpenlng  and  Pomeroy  shows  that  indeed  earth- 
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quake/exploslon  discrimination  can  be  performed  with  SH  wave  magnitudes. 
This  work  is  being  done  with  WWSSS  and  the  Canadian  network.  Such 
discrimination  should  be  attempted  with  long  period  LASA  data. 

However,  the  SH  wave  work  being  done  now  by  Turpening  and  Pomeroy 
is  restricted  to  large  explosions  (Benham,  Faultless,  Boxcar,  and 
Greeley)  and  comparable  sized  earthquakes.  The  problem  lies  with  the 
explosions.  Since  they  generate  relatively  little  SH  wave  radiation 
it  is  very  difficult  to  find  stations  situated  at  distances  large 
enough  to  separate  the  "weak"  SH  waves  from  the  surface  wave  train 
yet  operating  at  magnifications  large  enough  to  provide  reasonable 
measurements.  Because  of  this  problem  the  work  of  Turpening  and 
Pomeroy  cannot  be  applied  to  smaller  magnitude  explosions. 

Therefore,  it  is  recomnended  that 

1.  Long  period  LASA  data  be  used  to  study  the  discrimination 
capabilities  of  SH  wave  magnitude  vs.  P  wave  magnitude  diagrams. 

2.  Using  the  coiTq)uter/LASA  playback  system  developed  at 
the  Geophysics  Laboratory,  investigate  array  and/or  three-component 
signal  processing  techniques  to  enhance  weak  SH  arrivals  from  small 
explosions.  Such  a  study  should  make  use  of  not  only  LASA  data  but 
also  NOSAR  and  ALPA  data.  ALFA  in  fact  lies  at  a  good  distance  from 
NTS  such  that  the  SH  arrivals  would  be  separated  from  the  surface  wave 
t  rains . 

3.4.  SPECIAL  STTIDIES 

3.4.1.  MICROEARTHQUAKE  INVESTIGATIONS.  The  local  seismic  activity 
in  three  active  fault  areas  some  60  to  100  km  from  the  JORUM  event 
showed  no  significant  difference  before  or  after  the  event,  thus  de¬ 
monstrating  that  the  shot  did  not  generate  any  earthquakes  in  the  area 
of  these  temporary  stations.  High  record  gains  were  used  and  the  sta¬ 
tions  were  operated  continuously  so  no  significant  events  should  have 
gone  undetected. 
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l.ocally  very  small  events  were  detected  at  all  stations  that 
hail  spectral  peaks  as  high  as  80  II2.  The  Groom  Lake  station  had 
the  highest  number  of  these  events.  At  this  site  the  occurrence 
of  these  events  showed  a  strong  diurnal  effect. 

3. A. 2.  DLGOin^LING  EXPIilRlMENT.  The  seismic  measurements  of 
niAMONO  DUST  indicated  that  considerable  decoupling  was  achieved. 

The  absolute  efficiency  of  the  decoupling,  however,  cannot  be  de¬ 
termined  with  any  confidence  without  additional  tests.  The  five 
stations  used  to  record  this  event  should  be  reoccupled  for  any 
additional  shots  In  this  area.  The  Ideal  case  would  be  to  fire  a 
tamped  nuclear  shot  of  the  same  yield.  This  would  allow  highly 
accurate  decoupling  measurements. 

Smaller  yield  tamped  HE  shots  and  a  nuclear  shot  of  the  same 
yield  fired  in  the  same  sized  cavity  without  the  heat  sink  would 
provide  very  useful  data  for  computing  the  efficiency  of  the  DIAMOND 
DUST  experiment. 

Current  studies  are  being  made  of  data  obtained  from  the 
U.S.C.  &  G.  S.  on  the  HUDSON  MOON  and  DIESEL  TRAIN  shots.  While 
these  are  larger  events,  they  were  recorded  at  comparable  distances 
to  our  DIAMOND  DUST  recording  stations  and  the  event  was  fired  in 
tuff.  Hence  the  data  will  be  useful  for  scaling  purposes  and  perhaps 
a  better  estimate  can  be  made  on  the  decoupling  efficiency. 

3. A. 3.  OTHER  STUDIES,  The  conclusions  on  the  other  studies  may 
be  found  in  the  published  reports  referenced  in  Appendix  I-III. 

A 

EQUIPMENT  PURCHASES  AND  DEVELOPMENT 

Listed  below  are  all  of  the  capital  items  of  equipment  purchased 
under  this  contract. 
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I  Sprengnether  Recorder  (Purchased  August  1%9) 

Equipment  developed  under  this  contr'cc  include  the  following: 
Intermediate-depth  well  seismic  monitoring  system 
LASA  data  converter 
Optical  data  analyzer 
Transistorized  seismic  preamplifiers 
DC  Line  amplifiers 
Polyfilters. 

The  initial  design  on  the  above  equipment  was  done  under  sponsor¬ 
ship  of  previous  AFOSR  contracts  but  further  developments,  improvements 
and  some  construction  was  accomplished  under  this  contract. 
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I.eou  Reiter,  "An  Investisation  Into  the  Time  Term  Method  in  Refraction 
Seismolo,i;v Ikill.  Seism.  Soc.  Am..  Vol.  60,  No,  1,  pp.  1-13,  February 
1‘170. 


The  time  term  method  of  seismic  refraction  data  analysis 
allows  an  areal  arrangement  of  shot  points  and  stations  and 
yields  directly  a  three-dimensional  representation  of  under¬ 
lying  structure.  This  method  was  Investigated  analytically 
and  through  numerical  experiments  to  determine  implicit 
sources  of  error  and  to  provide  guidelines  for  Its  proper 
use.  Error  terms  for  computed  refractor  depths  were  derived 
for  two  site  arrays  over  two  structures,  a  simply  dipping 
refractor  surface  and  a  symmetric  anticlinal  refractor  sur¬ 
face.  These  error  terms  proved  to  be  qualitatively  reliable 
guides  to  tho  performance  of  the  time  term  method  using 
models  of  multi  site  arrays  over  planar  and  blplanar  refrac¬ 
tor  structures.  Neither  continuity  of  time  terms  nor  use  of 
a  known  refractor  velocity  In  the  computations  offered  a 
guarantee  of  a  more  accurate  refractor  reconstruction.  In 
the  two  models  tested  a  refractor  velocity  greater  than  the 
true  velocity  yielded  the  best  reconstruction.  Sensitivity 
to  angle  of  dip  indicated  preferred  orientation  of  shot- 
station  lines  along  gentle  apparent  dips  when  the  underlying 
structure  is  known  and  their  orientation  along  several  azi- 
muths  when  the  structure  Is  unknown  so  as  to  mute  the  effect 
of  Iilgh  dip  angles.  If  the  structure  has  refractor  surfaces 
dipping  more  than  a  few  degrees,  the  time  term  method  will 
give  Inaccurate  results  unless  previous  knowledge  of  the  re¬ 
fractor  configuration  permits  model  studies  to  determine 
error  terms  and  arrangement  of  shots  and  stations.  Unfor¬ 
tunately  such  knowledge  Is  usually  the  goal  of  a  seismic 
experiment . 


C.  G.  bufe  and  D.  E.  Willis,  "A  Note  on  Short  Period  Seismic  Spectral 
Ratios,"  Earthquake  Notes.  Vol,  40,  No,  1,  March  1969, 

An  inverse  square  frequency  dependence  of  the  amplitude 
ratio  of  P  waves  from  events  differing  in  yield  by  several 
orders  of  magnitude  is  shown  to  exist  over  several  octaves 
in  frequency.  The  behavior  of  the  spectral  ratio  of  parti¬ 
cle  velocities  for  the  12.5  kiloton  SHOAL  nuclear  explosion 
and  2.5  ton  USGS  Delta,  Utah  shot  is  used  to  estimate  the 
radius  of  "equivalent  radiator"  at  SHOAL. 
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P,  L,  Jackson,  "Seismic  Ray  Simulation  for  Spherical  Earth,"  Bull,  Seism. 

Soc,  Am.  .  Vol,  60,  No,  3,  pp,  1021-1025,  June  1970. 

This  letter  describes  a  method  of  ray  simulation  to  relate 
velocity  distribution  to  travel  times  in  spherical  earth  which 
also  is  based  upon  Snell's  law  of  refraction  with  the  addition 
of  the  law  of  reflection.  Given  a  velocity  distribution,  tra¬ 
vel  times  for  all  exclusively  P  or  S  phases  of  body  waves  can 
be  generated  and  a  predetermined  number  of  multiple  reflections 
included.  No  restrictions  exist  upon  the  positive  or  negative 
values  of  dv/dr  or  upon  discontinuities  in  the  velocity  dis¬ 
tribution.  Tables,  mathematical  functions,  or  combinations  of 
both  are  used  to  refer  to  the  velocity  distribution.  Mathe¬ 
matical  functions  can  be  used  within  depth  ranges  where  it  is 
convenient  to  represent  the  distribution  mathematically;  tables 
in  depth  ranges  with  abrupt  or  anomalous  distributions,  or 
where  one  desires  to  alter  specific  velocity  values  when  re¬ 
peating  the  tracing  of  a  ray. 


D.  E.  Willis  and  J.  T.  Wilson,  "A  Note  on  the  Anna,  Ohio  Earthquake  of 
July  26,  1968,"  submitted  for  publication  to  Earthquake  Notes.  July  1970, 

A  small  earthquake  located  near  Anna,  Ohio  was  recorded 
on  the  University  of  Michigan's  special  well  system  on  July 
26,  1968.  The  body  wave  magnitude  was  determined  to  be  3.0 
with  a  modified  Mercalll  intensity  of  at  least  III.  Spectral 
analyses  showed  that  most  of  the  seismic  energy  was  contained 
in  a  band  from  3  to  10  Hz. 
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Appendix  II 

CllUONOLOcM  C/\I,  ANNOTATED  BIBLIOGRAPHY  OF  TECHNICAL  REPORTS 
SPONSORED  bY  AFA9(638)-1759 


David  E.  Willis  and  Philip  L.  Jackson,  Collection  and  Analysis  of  Seismic 
Wave  Propagation  Data  (Annual  Report),  Rept.  No.  8071-6-P,  Willow  Run 
Laboratories  of  the  Institute  of  Science  and  Technology,  Geophysics 
Laboratory,  The  University  of  Michigan,  Ann  Arbor,  August  1967. 

This  report  summarizes  one  year  of  theoretical  and  applied 
research  on  propagation  of  seismic  waves  and  techniques  for 
analyzing  data.  The  main  objectives  were  to  determine  the  fre¬ 
quency  and  energy  of  seismic  signatures,  and  investigate  atten¬ 
uation,  patterns  of  azimuthal  radiation  from  source  regions, 
and  methods  of  determining  the  type  of  motion  at  the  source. 

Natural  and  artificial  sources  were  studied  to  develop  diagnos¬ 
tic  aids  for  distinguishing  between  earthquakes  and  underground 
nuclear  detonations.  Equipment  for  selection,  reformatting, 
and  dlgltal-to-analog  conversion  for  digitally  recorded  LASA 
data  was  constructed  and  is  being  checked  out.  Several 
approaches  for  using  the  parallel  computational  capabilities  of 
optics  for  LASA  data  were  developed.  A  study  of  background 
noise  and  reciprocity  for  teleseismic  events  as  recorded  on  the 
bottom  of  a  large  fresh  water  lake  has  commenced  with  the  em¬ 
placement  of  three- component  seismometers  in  Lake  Superior* 

Array  data  have  been  used  for  crustal  studies  on  the  Eastern 
United  States.  Digital  mode  filtering  was  investigated*  A 
perturbation  theory  for  seismic  sources  was  developed* 


David  E*  Willis,  An  Investigation  of  Seismic  Wave  Pronaeation  in  the 
Eastern  United  States.  Rept*  No*  8071-16-T,  Willow  Run  Laboratories 
of  the  Institute  of  Science  and  Technology,  Geophyaics  Laboratory, 
The  University  of  Michigan,  Ann  Arbor,  July  1968* 

This  paper  describes  the  travel-time  anomalies  and  atten¬ 
uation  losses  of  seismic  compressional  waves  generated  by  a 
series  of  underwater  explosions  in  the  Eastern  United  States* 

The  efficient  tamping  of  the  shots  fired  in  water  provided  a 
seismic  source  that  could  be  detected  at  much  larger  ranges 
than  could  he  accomplished  by  equivalent  sized  shots  fired 
underground* 

A  number  of  mobile  field  recording  stations  equippad  with 
three-component  matched  short  period  seismometers  and  ZMgnatic 
tape  recorders  were  used  to  record  273  individual  shots.  A 
total  of  1,293  recordings  tiere  obtained  along  a  rsversad  pro¬ 
file  extending  from  International  Falls,  Mlnnasota  to  the 
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Atlantic  Coast.  The  underwater  shots  were  fired  In  Lake 
Superior  and  the  Atlantic  Ocean.  Precise  travel  times  were 
obtained  by  recording  radio  time  signals  at  all  of  the  re¬ 
cording  stations. 

An  analysis  of  the  travel-times  of  the  seismic  waves 
disclosed  that  the  earth’s  crust  varies  In  thickness  from 
28.1  km  near  the  Atlantic  Coast  in  North  Carolina  to  50.4 
km  near  the  Keweenaw  Peninsula  in  upper  Michigan.  The  crust 
in  North  Carolina  was  found  to  be  comprised  of  a  single  lay¬ 
er  with  a  compressional  wave  velocity  of  6.0  km/sec.  A  two- 
layer  crust  with  compressional  wave  velocities  of  6,2  and 
7.7  km/sec  was  disclosed  in  upper  Michigan.  Travel  time  re¬ 
siduals  across  the  Appalachian  Mts  would  indicate  a  mountain 
root  system  similar  to  that  found  under  the  Rocky  Mts. 

The  Lake  Superior  area  was  found  to  be  more  efficient  in 
the  coupling  of  energy  into  seismic  waves  by  the  underwater 
shots  than  in  the  Atlantic  Ocean.  This  is  believed  due 
mainly  to  the  bubble  pulse  phenomenon  and  the  reinforcement 
caused  by  reverberation  of  sound  waves  between  the  lake 
bottom  and  the  water's  surface. 

A  new  technique  of  displaying  seismic  attenuation  data  in 
three  dimensional  space  was  developed  that  permits  presen¬ 
tation  of  large  quantities  of  data  in  a  concise  manner,  hence 
aiding  in  the  interpretation  of  the  data.  Consistent  anoma¬ 
lies  in  the  spectrum  of  the  seismic  waves  were  f  und  that 
could  be  correlated  with  various  parameters  in  tne  source 
region  and  with  propagation  effects. 

The  detailed  spectral  studies  that  were  made  as  a  function 
of  distance  from  the  source  disclosed  that  the  attenuation  of 
the  seismic  waves  could  be  expressed  by  the  equation 


where  A  ■  amplitude,  X  -  distance,  n  >  geometrical  spreading 
factor,  a  ■  absorption  coefficient  and  f  ■  frequency.  The 
crustal  attenuation  data  in  the  Lake  Superior  area  were  found 
to  have  a  value  of  n  ■  1.5  and  a  ■  2-3  x  10“3,  Upper  mantle 
compressional  wave  attenuation  data  gave  a  value  of  n  ■  1  and 
a  -  2  X  10-^, 


David  E.  Willis  and  Philip  L.  Jackson,  Collection  and  Analysis  of  Seismic 
Wave  Propagation  Data  (Annual  Report),  Rept.  No.  8071-15-P,  Willow  Run 
Laboratories  of  the  Institute  of  Science  and  Technology,  Geophysics 
Laboratory,  The  University  of  Michigan,  Ann  Arbor,  September  1968. 

This  report  discusses  the  results  of  an  extensive  long- 
range  reversed  refraction  profile  which  traverses  the  Michigan 
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Basin  and  the  Appalachian  Mountains.  Particular  emphasis 
is  placed  on  the  attenuation  of  first  compressional  wave 
arrivals  and  on  crustal-upper  mantle  structure.  Also  in¬ 
cluded  are  results  of  a  survey  on  lake  bottom  seismic  back¬ 
ground  noise  in  Lake  Superior.  Although  large  signal  levels 
were  obtained  in  the  lake  bottom,  long  periods  of  high  back¬ 
ground  noise  indicated  that  land  based  seismographs  are 
superior.  Theoretical  studies  in  elastic  wave  scattering 
were  made  to  evaluate  the  effect  of  nonhomogeneities  on  the 
propagation  of  elastic  waves.  Theoretical  development  of 
mode  filtering,  (previously  demonstrated  empirically),  is 
nearly  completed.  A  LASA  digital  data  convertor  was  de¬ 
signed,  constructed,  and  put  into  operation.  Two-  and  three- 
dimensional  seismic  ray-tracing  techniques  using  digital 
modeling  were  developed.  These  techniques  enable  the  seis¬ 
mologist  to  determine  the  effects  of  lateral  Inhomogeneltles 
and  irregularities  in  velocity  Interface  on  travel  time  and 
attenuation  of  seismic  waves. 


Philip  L.  Jackson,  Roger  M.  Turpenlng,  and  David  E.  Willis,  Collection 
and  Analysis  of  Seismic  Wave  Propagation  Data  (Annual  Report),  Kept. 
No.  a071-25-P,  Willow  Run  Laboratories  of  the  Institute  of  Science  and 
Technology,  Geophysics  Laboratory,  The  University  of  Michigan,  Ann 
Arbor,  July  1969. 

The  work  performed  under  this  contract  during  the  past  year 
Included  mode  filtering,  S-wave  studies;  seismic  ray  tracing; 
theoretical  studies  in  elastic-wave  scattering;  seismic  wave 
propagation  near  the  core-mantle  boundary;  time  term,  travel  time, 
spectral,  and  attenuation  studies;  and  field  neasurements  pri¬ 
marily  with  an  intermediate-depth  well  seismograph  station. 

Seven  Journal  articles  were  either  published  or  submitted  for 
publication  during  this  period.  Continuation  of  this  contract 
will  be  concerned  with  S-wave  studies  for  difference  in  arrival 
times  caused  by  layering  at  the  recording  location,  simulation 
by  ray  tracing  to  investigate  the  velocity  structure  underlying 
LASA,  and  a  field  measurement  program  of  continuous  recordings 
in  Nevada. 


David  E.  Willis,  Philip  L.  Jackson,  and  Roger  M.  Turpenlng,  Collection 
and  Analysis  of  Seismic  Wave  Propagation  Data  (Semi-Annual  Report) , 

Rept.  No.  8n71-3n-L,  Willow  Run  Laboratories  of  the  Institute  of  Science 
and  Technology,  Geophysics  Laboratory,  The  University  of  Michigan,  Ann 
Arbor,  January  1970. 
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The  scope  of  this  contract  Includes  S-wave  mode  filtering, 
seismic  ray  tracing,  and  microearthquake  investigations.  After 
much  delay  pertinent  LASA  tapes  have  been  obtained;  the  LASA 
data  processor  has  been  interfaced  with  the  PDP-8/I  computer 
located  in  the  Geophysics  Laboratory  to  perform  the  S-wave 
mode  filtering. 

Highly  accurate  seismic  ray  simulation  through  two-dimen¬ 
sional  heterogeneous  velocity  distributions  has  been  achieved, 
and  the  extension  to  three  dimensions  with  selective  regions 
for  sampled  or  mathematically  represented  data  is  being  pro¬ 
grammed,  Preliminary  analysis  of  seismograms  recorded  before 
and  after  the  JORUM  underground  nuclear  explosion  indicate  that 
there  was  no  significant  difference  in  the  occurrence  of  local 
microearthquakes  before  and  after  this  explosion. 


Charles  G.  Bufe,  An  Estimate  of  the  Configuration  of  the  Surface  of  the 
Earth *8  Core  from  the  Consideration  of  Surface  Focus  PeP  Travel  Times. 
(PhD  Dissertation)  and  (Supplement  to  8071  Final  Report),  Rept.  No. 
8071-33-F2,  Willow  Run  Laboratories  of  the  Institute  of  Science  and 
Technology,  Geophysics  Laboratory,  The  University  of  Michigan,  Ann  Arbor, 
July  1970. 

Travel  times  of  PeP  and  P  phases  from  nuclear  and  high 
explosive  sources  are  interpreted  in  terms  of  variations  in 
the  radius  of  the  earth's  outer  core.  Travel-time  residuals 
from  the  Taggart-Engdahl  (1968)  tables  are  determined  for  the 
PeP-P  time  interval  and  for  PeP  travel  times,  corrected  for 
station  and  source  terras.  The  travel  tiraes  are  also  corrected 
for  elevation  of  source  and  station  and  ellipticity,  excluding 
the  core  ellipticity  terra.  This  interpretation  favors  a  core 
which  is  slightly  larger  than  that  of  the  reference  raodel  of 
Taggart  and  Engdahl  (1968)  and  which  has  less  ellipticity  than 
estlraated  by  Bullen  (1936).  Other  interpretations  in  terras  of 
lateral  variation  in  mantle  velocity  are  also  possible. 

Because  all  of  the  events  and  most  of  the  stations  whose 
data  are  used  lie  in  the  northern  hemisphere,  the  description 
of  core  shape  is  Incomplete  and  is  representative  only  of  re¬ 
gions  of  the  core  surface  in  the  vicinity  of  the  PeP  reflection 
points. 

The  distribution  of  data  does  not  permit  a  conventional 
spherical  harmonic  analysis  of  the  Inferred  variation  in  core 
radius.  A  modified  spherical  harmonic  analysis  method  is  de¬ 
vised  to  smooth  the  data  and  to  estimate  the  shape  of  the  core. 

The  variation  of  core  radius  determined  from  this  representation 
is  approximately  i:10  kilometers  on  the  basis  of  terms  to  degree 


19 


INMIII  II  tn  SI  UNI  I  AM)  III  IINOl.lHiV 
MM  I  Nl\  I  Ksn  A  1)1  MU  nil. AN 

anil  order  ').  The  standard  deviation  of  the  data  is  approxi¬ 
mately  *20  kilometers,  indicating  the  presence  of  variations 
of  higher  degree  and  order,  many  of  which  are  probably  due  to 
factors  otlier  tlian  undulation  of  the  mantle-core  boundary. 
Variations  in  geoid  height  resulting  from  the  inferred  vari¬ 
ations  in  Cl' re  radius  do  not  correlate  with,  but  are  of  the 
same  order  of  magnitude  as,  those  determined  from  satellite 
observations . 

/\nalysis  of  a  limited  number  of  short-period  PcP  amplitudes 
shows  no  significant  variation  of  the  PcP  amplitude/period 
ratio  with  eplcentral  distance  and  an  estimated  scatter  of  data 
on  tlie  same  order  as  that  found  for  P  in  the  distance  range  70® 
to  90®.  The  PcP  phase  is  generally  shorter  in  apparent  period 
than  tlie  P  phase.  This  relationship  is  reversed,  however,  for 
arrivals  from  the  Novaya  Zemlya  event  of  October  27,  1966, 

The  contribution  of  a  hypothetical  upper  mantle  low  velocity 
zone  of  varying  thickness  to  the  amplitude  and  travel-time 
scatter  is  demonstrated. 


Philip  L.  Jackson,  Digital  Simulation  of  Seismic  Rays  (PhD  Dissertation) 
and  (Supplement  to  8071  Final  Report),  Rept.  No,  8071-33-F2^,  Willow  Run 
Laboratories  of  the  Institute  of  Science  and  Technology,  Geophysics 
Laboratory,  The  University  of  Michigan,  Ann  Arbor,  July  1970, 

Simulation  of  seismic  rays  for  a  spherical  earth  and  a  flat 
earth  has  been  achieved  in  highly  complex  models.  Travel  times 
and  approximate  amplitudes  of  seismic  waves  can  be  found  for 
both  two-  and  three-dimensional  models  of  portions  of  the  earth. 
Hypotheses  of  complicated  and  detailed  velocity  distributions' 
within  the  earth  can  be  tested  by  comparing  travel  times  found 
in  models  with  travel  times  found  from  actual  earthquakes  and 
explosions. 

Since  ancient  times  rays  have  been  constructed  to  denote  the 
direction  in  which  energy  is  propagating.  As  in  other  scienti¬ 
fic  fields,  the  construction  of  rays  to  represent  propagation 
paths  has  been  a  fundamental  technique  of  modem  seismology  since 
its  inception  in  the  latter  part  of  the  19th  century.  In  seis¬ 
mology  and  other  disciplines  ray  construction  has  been  applied  to 
simplified  geometries.  It  has  been  necessary  to  assume  that  the 
seismic  wave  velocity  distribution  of  the  earth  was  relatively 
uniform  and  symmetric. 

Recently,  however,  the  earth  has  been  found  to  be  more  complex 
and  non-uniform  than  formerly  assumed.  Travel  times  have  been 
extensively  measured  with  more  accurate  instrument  techniques, 
making  historical  methods  of  ray  construction  inadequate  for 
either  visualization  of  the  rays  or  numerical  indication  of 
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expected  travel  times,  A  need  has  thus  arisen  in  seismology 
to  test  highly  heterogeneous  models  of  seismic  velocity  distri¬ 
bution,  At  the  same  time  the  development  of  the  modern  digital 
computer  has  provided  a  means  of  performing  the  necessary  ray 
constructions  and  numerical  calculations. 

The  problem  of  complicated  seismic  velocity  distributions 
was  therefore  Investigated  in  terms  of  the  most  appropriate  use 
of  the  digital  computer.  For  this  Investigation  a  velocity 
field  was  set  up,  and  the  propagation  computations  made  for 
short  segments  of  rays  within  this  field.  For  each  segment  the 
surrounding  parameters  of  propagation  were  known,  hence  changes 
in  direction  and  time  of  travel  could  be  computed,  and  newly 
computed  segments  joined  to  those  preceding  until  a  complete 
ray  is  constructed.  The  total  travel  times  are  also  found  by 
adding  the  travel  times  of  all  the  joined  ray  segments.  Essen¬ 
tially,  the  nature  of  propagation  was  duplicated  on  the  computer, 
in  that,  at  the  location  of  each  segment  along  the  path  of  pro¬ 
pagation,  the  initial  condition  and  effect  of  the  surroundings 
determine  the  succeeding  direction  of  the  following  segment. 

Both  visual  and  numerical  results  have  shown  that  this  simu¬ 
lation  method  can  be  usefully  applied  to  investigation  of  seismic 
velocity  distributions  of  portions  of  the  earth  of  any  size  or 
complexity.  Computer  drawn  plots  of  the  entire  spherical  earth, 
and  on  vertical  cross-sections  of  the  underlying  structure  across 
California  and  of  typical  lithospheric  plate  boundaries  Illustrate 
useful  applications  of  this  method.  Accuracy  is  illustrated  by 
calculation  of  essentially  the  same  travel  times  for  P  waves  as 
published  by  the  Selsmologlcal  Society  of  America, 


David  E,  Willis  and  Robert  F,  Hand,  Seismic  Measurements  of  DIAMOND  DUST 
(Supplement  to  8071  Final  Report),  Report  No,  8071-33-F3/3620-1-T, 

Willow  Run  Laboratories  of  the  Institute  of  Science  and  Technology,  The 
University  of  Michigan,  Ann  Arbor,  July  1970. 

A  small  nuclear  device  detonated  in  an  underground  cavity 
was  recorded  at  five  close-in  sites.  Portable  stations  were 
located  to  study  azimuthal  effects.  Strong  asymmetries  were 
observed  in  frequency  content,  amplitude  and  wave  train  dura¬ 
tion,  These  effects  can  be  explained  in  part  by  the  variations 
in  geology  at  the  recording  sites.  However,  source  radiation 
could  account  for  some  of  the  observed  anomalies. 
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n.  E.  Willis  and  J.  T.  Wilson,  "A  Long-Range  Reversed  Seismic  Profile  in 
the  Eastern  United  States,"  presented  at  the  International  Union  of 
Ceodesy  and  Ceophysics,  XIVth  General  Assembly,  Switzerland,  September 
25  through  October  7,  19b7. 

A  series  of  underwater  high-explosive  shots  fired  during 
the  past  five  years  have  provided  data  for  a  reversed  seismic 
profile  approximately  2000  km  in  length.  These  shots  include 
the  19b2  offshore  North  Carolina  experiment,  the  1965  ECOOE 
series,  the  CHASE  III,  IV,  and  VII  events,  and  the  Lake 
Superior  shots  of  1963,  1964,  and  1966,  Close  station  inter¬ 
vals  were  obtained  along  the  profile.  Travel-time  curves, 
frequency  analyses  and  attenuation  measurements  are  presented 
for  these  recordings.  An  intermediate  crustal  layer  with  a 
velocity  of  7.7  km/sec  was  disclosed  on  the  northwestern  por¬ 
tion  of  this  profile.  Higher  Pn  velocities  and  a  thicker  crust 
was  also  indicated.  Smaller  travel  time  residuals  across  the 
Michigan  Basin  were  observed  on  both  profiles.  The  travel-tine 
data  indicated  a  possible  mountain  root  system  under  the 
Appalachians.  Contoured  particle-velocity  maps  made  fron  the 
spectral  data  disclosed  amplitude  anomalies  that  could  be  cor¬ 
related  with  conditions  in  the  source  region.  These  results 
aid  in  the  interpolation  of  the  attenuation  data. 


C,  G.  Bufe  and  D.  E.  Willis,  "Seismic  Radiation  from  SHOAL,"  presented 
at  the  64th  Annual  Meeting  of  the  Selsmological  Society  of  America, 
Tuscon,  Arizona,  April  11-13,  1968, 

Travel  times  and  particle  velocities  have  been  examined  for 
distances  to  1600  km  and  for  widely  distributed  azimuths  for 
the  SHOAL  event.  A  composite  curve  of  raw  travel-time  data  re- 
affirms  the  anomalously  low  P^  velocities  previously  observed 
in  the  Basin  and  Range  province  and  yields  a  crustal  thickness 
of  about  26  km.  Particle  velocities  were  measured  for  3  sets 
of  )  stations  each  at  distances  of  250,  320,  and  approximately 
1100  km.  The  particle  velocity  spectre  vary  widely  from  site 
to  site  and  generally  refl'^ct  characteristics  of  the  spectra 
of  the  seismic  background  nr  e^ich  site. 
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David  K.  UllllH,  "Selamlc  Nolac  on  ch«  Bottom  of  I.ak«  Superior «"  proKented 
at  the  49th  Annual  Meeting  of  the  American  Ceophyeical  Dnion,  Uanhlngton, 
1).  r, ,  April  8-11,  1968. 

A  RcrloR  of  Reiemic  hackgroiind  nciMw  meaauremente  were  made 
on  the  bottom  of  Lake  Superior  diirintt  .iune-AuRuut  1967.  The 
mcaRurementfl  were  made  at  the  eite  of  the  1966  lake  Superior 
Rhut  point  In  approximately  S75  feet  of  water.  Two  inetrunent 
packages  that  were  obtained  on  loan  from  the  I’hivoniity  of 
California  ware  ueed  for  theae  meaaurementa*  Lach  was  equip¬ 
ped  with  three-component,  ahort-period  aeie*%imeterN,  a  magnetic 
tape  recorder,  and  en  accurate  chronometet.  (Jne  free-fall  drop 
and  throe  tethered  dropa  were  made.  Approximately  ten  deye  of 
contiauoua  recordlnge  were  obtained  for  each  of  three  time 
perloda,  the  laat  of  tfhlch  Included  a  tethered  unit  and  a  free- 
fall  unit.  The  background  nolae  wan  found  to  be  quite  high, 
averaging  In  the  vicinity  of  100  mu  at  1  cpa.  Long  perloda 
loatlng  up  to  aeveral  dayo  of  extremely  high  background  nolae 
were  noted.  Theae  perloda  could  be  correlated  with  wave 
action  and  meteorological  conditions*  At  their  peaka,  ground 
dloplacemanto  of  the  order  of  1000  mu  were  obaerved.  Several 
earthquakea  were  reoordad  during  thla  period  ae  well  aa  many 
quarry  ahota.  Spectral  analyoea  of  tliMe  eventa  are  coirrelated 
with  data  obtained  from  a  etatlon  operated  nearby  on  land. 

Coiparlaona  of  algnal-to-nolae  ration  are  alao  preaented. 


Philip  L.  Jackaon,  "TVo-  and  Threo-tllmeoalonal  Ray  Tracing,  Travel  fine, 
and  Aaplltuda  Conpucaklon,**  preaented  at  the  hith  Annual  Meeting  of  the 
Selamologlcal  Society  of  America,  Tuacon,  Arlions,  April  li-*l,  I9M, 

Selomlc  rryn  can  be  traced  through  velocity  end  denalty 
atructurea  of  any  degree  of  Inhomagemalty.  A  digital  co» 
puter  la  uned  to  trace  and  graphically  to  preaent  the  raya, 
and  to  compute  the  travel  tinea  end  anplltudea  of  theae 
emerging  at  the  aurface.  The  mother  of  multiple  reflectloma 
can  bn  aelected  or  con  be  limited  by  a  mlnlmnn  amplitude  upon 
reflection.  A  aompllng  of  the  atmeture  on  two-  or  three- 
dlnanelenal  grlda  eenprlaea  the  data  for  the  ray  treeing.  An 
the  data  can  be  "pat^d**  tofither  and  the  cemtinuatlon  of 
raya  traced  rough  each  data  "patch,"  only  a  amall  ameuni  of 
cenputer  memory  la  requlredi  thma,  detail  la  not  llmitad  by 
eenputer  mamory.  Caleonp  plot*  of  ray  tracea  through  aeveral 
repreaeotative  velocity  atructurea  have  been  made.  Inclmled 
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Iwilat  ta  la  l»M  la  Wll»iiia»  taaiai at  ll•IA•  lam* 

A  trlAartlia  aat  af  tlini  lanpaiyiat  paHaAla  tlaftfarlaA 
aaliMtf^  AtatlaM  mat  aptfatai  faatlaaaaily  far  ifna  aaaia 
ptlar  aM  lAnaa  wtata  aflat  ila  JMIl  MaltftfMMl  aatlaar 
iiat  af  iffiaiAtf  ll«  lAM*  IlMaa  atailaaa  aaia  latai«i  M 
la  IM  ia  ffM  Kka  mmm  aaat  attlaa  faali**  Ma  al^ilfl- 
«Mit  aiMiiiaa  mm  iAatmA  la  laari  Mctiaalaalt  a«.ilaltf 
iaftata  mi  aflat  tl»  aami  laiicailaa  iIm  tim  Al4  aai  cawia 
aay  talaaM  la  latiaalc  aatm  la  iila  alclaliy  af  tlw  taaparary 
•laUiM*  Hw  wIdiiaN  ftaawaty  taifMa  af  tia  tacardlai 
aftiami  allawN!  ilta  dMacilaa  af  fatal  ataaia  lAai  caaialaad 
ai#  fmiaatf  aaatat  ailtA  paiAiad  aa  Algl  aa  it*  Ai  aaa 
aiailaa  a  ptaaaaatad  dlataal  affaci  la  lia  atcartiaca  af 
uwaa  hi#  ftapaMMf  aaaaia  aaa  ahitrtal. 


•avid  fc«  milla*  "Aitamhal  Affwatrlaa  flhaaiaai  Araaad  a  halat  ioarca/ 
fat  ptaaiiailaa  ai  iht  diAfSiA  Aaaaal  Umlaa  ta  ha  hald  la  Hi  laaidtat» 
iaaMdat  II•1A•  If  Id, 

lha  daafhfalca  Uhataiary  ai  lha  tlilvaraliy  af  Hlciil«aa 
apataiad  fiva  aiatt>farlad  aalaaatt#!  atailaaa  la  racatd  a 
laa-ylald  aadtttiaaad  aatlaat  ahai  fitad  at  tha  havada  that 
ilia,  IW  ahai  aai  fitad  la  m  alt^flllad  cavity  la  taat 
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till'  1  iMs  ill  i  1  i  tv  of  tliemio-nioc-lianlcal  tei.'lmiqui>K  for  de- 
ooiipliiu*,  siismii-  fiUMitv,  Four  stations  wore  located  at 
.ipt>  I  oxiiiiati'l  V  kn  to  study  the  aititnutlial  effect,  Tlio 
tiftli  station  was  located  at  14  km  to  provide  attenua¬ 
tion  data.  Strong  asymmetries  were  oliserved  in  frequency 
(-•ontent,  amplitude  and  wave  train  duration.  Tlie  complicn- 
teil  local  v.eoli-’Ky  could  explain,  in  part,  the  ohsorved 
anomalies.  However,  at  the  present  time,  some  source  radi¬ 
ation  effects  cannot  be  discounted. 
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A  anr  mcUnd  of  aioMlailait  otlOMle  nya,  eonpaiiag  travtl  tlnw,  ud  ■uproalatlat  aapUladM  la  aaita 
Mdtla  baa  bttn  dowloiwd.  Tba  diftlal  conpaitr  Is  aaod  to  aodal  a  apbarteal  or  a  flat  aaitb.  1>ro- or  tbroo- 
dlawaa tonal  modrla  art  tmpktytd.  Tha  oNthod  can  bt  nood  for  batarofiBaona  nwdala,  and  laeladaa  cowpniailor. 
of  awblpla  raflactionB.  Traaal  unwa  bata  baan  obtalaad  for  a  nwdal  of  a  craan  aactlon  ondarlylat  LAIA. 

la  INT  tha  lavastiftaiion  of  optical  procasalat  oaa  taradaBtad.  bdtlaiad  la  Ml,  tba  Inaaotlsattoa  oaa  coa> 
dactad  to  aaplora  tha  aMicatbia  of  eobarani  optical  procaaolag  to  aatsaUe  aaalpols.  Da  aaa  oaa  foaad  to  lla 
prtauirlly  In  oil  aiploratlan  salamologjr. 

A  data  proeasstnii  capabiiay  apaetflcaUy  lataadad  for  tba  stady.of  Iba  ••oaaa  paitlela  aarlattoaa  aeroos 
LAIA  oaa  bvilt  on  ibla  contract.  Tba  LAIA  dau  proeasaor,  obicb  aalstad  la  tbia  laboratorr  for  tha  optical  I 
aaalysaa  of  LAIA  data,  was  aaNllflad  to  latarfaca  with  this  laboratory’s  PDP*I  coavolar.  This  data  procasalnp 
capability  was  n  t  uaad  dot  to  foadlnf  IlmltatlonB. 

lOcrosatamic  invamieatlona  In  conaactloa  with  JORUM  ladlcatad  that  tha  flrlnp  of  a  larfs  ytald  aadargroo 
ahot  did  riot  itroduca  any  mvaaurabla  tactonic  raiaosa  at  anarcy  la  actira  taalt  raflooa  at  a  distanca  of  60  to  100  | 
bm  from  the  iniirce.  S**|»mic  maasuramants  of  DIAMO^»D  DUIT  Indicated  dacoa^ap  was  achlavad  but  tha  ac> 
curate  dairrmination  of  iha  afficlancy  of  tha  dacoapUap  will  haaa  to  ba  dataraUaad  by  additional  taata.  A  pro- 
aounced  asymmatry  in  aciamlc  wave  charactarlstlca  was  obaanrad. 

Travel  time  studies  of  PeP  and  P  wart  used  to  datarmlna  tha  aorfaca  coofifaratlon  of  a  portion  of  tha 
earth's  core.  An  extensive  study  was  made  on  tha  attenuation  and  travel  IIrmb  of  salamic  waves  la  tha  Eastern 
United  States.  Theoretical  studies  wart  made  to  datarmlna  the  scattarlnp  affacta  of  plana  elastic  waves  by  sur¬ 
face  imperfect  ions  and  lo  determine  the  s  nircaa  of  error  la  tha  tima-tarm  method  In  refraction  salsmolo^. 

Spectral  studies  were  made  of  undersround  nuclear  shots  and  aarthquakas  recorded  at  talasalsmic,  ra- 
ploKit  and  near  distances.  An  additional  study  was  made  of  tha  charactarlstlca  of  seismic  backiprouiid  noise  re¬ 
corded  on  Ml'’  ifittiim  Ilf  Liiio  Sbnerior. 
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